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Superconductivity in a Two-Orbital Hubbard Model with Electron and Hole Fermi 
Pockets: Application in Iron Oxypnictide Superconductors 
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We investigate the electronic states of a one-dimensional two-orbital Hubbard model with 
band splitting by the exact diagonalization method. The Luttinger liquid parameter K p is 
calculated to obtain superconducting (SC) phase diagram as a function of on-site interactions: 
the intra- and inter-orbital Coulomb U and U' , the Hund coupling J, and the pair transfer 
J'. In this model, electron and hole Fermi pockets are produced when the Fermi level crosses 
both the upper and lower orbital bands. We find that the system shows two types of SC 
phases, the SC I for U > U' and the SC II for U < U' , in the wide parameter region including 
both weak and strong correlation regimes. Pairing correlation functions indicate that the most 
dominant pairing for the SC I (SC II) is the intersite (on-site) intraorbital spin-singlet with 
(without) sign reversal of the order parameters between two Fermi pockets. The result of the 
SC I is consistent with the sign-reversing s-wave pairing that has recently been proposed for 
iron oxypnictide superconductors. 

KEYWORDS: iron oxypnictide superconductors, two-orbital Hubbard model, pairing symmetry, exact diagonaliza- 
tion 



1. Introduction 

The recent discovery of iron oxypnictide superconduc- 
tors 1-5 with transition temperatures of up to T c ~ 55K 
has stimulated much interest in the relationship between 
the mechanism of the superconductivity and the orbital 
degrees of freedom. First-principles calculations have 
predicted the band structure with hole Fermi pockets 
around the T point and electron Fermi pockets around 
the M point. 6-8 By using the weak-coupling approaches 
based on multiorbital models, the spin-singlet s-wave 
pairing is predicted, where the order parameter of this 
pairing changes its sign between hole and electron Fermi 
pockets (sign-reversing s-wave pairing). 9-13 This uncon- 
ventional s-wave pairing is expected to emerge owing to 
the effect of antifcrromagnetic spin fluctuations. Since 
the strong correlation between electrons is considered to 
play an important role in the superconductivity of iron 
oxypnictides as well as in that of high-T c cuprates, non- 
perturbative and reliable approaches are required. 

As a nonperturbative approach, the exact diagonal- 
ization (ED) method has been extensively applied in the 
Hubbard, d-p, and t-J models. 14 ' 15 Although these mod- 
els are much simplified and mostly limited to one di- 
mension, it has elucidated some important effects of a 
strong correlation on superconductivity Using the ED 
method, we have studied the one-dimensional (ID) two- 
orbital Hubbard model in the presence of the band split- 
ting A. It is found that the superconducting (SC) phase 
appears in the vicinity of the partially polarized ferro- 
magnetism when the exchange (Hund's rule) coupling J 
is larger than its critical value on the order of A. 16 The 
result suggests that spin triplet pairing emerges owing 
to the effect of ferromagnetic spin fluctuation. In the 
case of A = 0, spin triplet superconductivity has also 



been discussed on the basis of bosonization 17-19 and nu- 
merical 20-22 approaches. Previous works, however, were 
restricted to the case of a single Fermi surface, and the 
effects of electron and hole Fermi pockets on supercon- 
ductivity have not been discussed therein. 

In this study, we investigate the ID two-orbital Hub- 
bard model with electron and hole Fermi pockets, where 
the Fermi level crosses both the upper and lower bands 
in the presence of a finite band splitting A. Using the 
ED method, the Luttinger liquid parameter K p is calcu- 
lated to obtain the SC phase diagram as a function of 
on-site Coulomb interactions in a wide parameter region 
including both weak- and strong-correlation regimes. It 
would clarify the effects of a strong correlation on su- 
perconductivity in iron oxypnictides. We also calculate 
various pairing correlation functions and discuss a pos- 
sible pairing symmetry. Although our model is much 
simplified and limited to one dimension, we expect that 
the essence of the superconducting mechanism of iron 
oxypnictides can be discussed. 

2. Model and Formulation 

We consider the one-dimensional two-orbital Hubbard 
model given by the following Hamiltonian: 



H = t ^ (4,m,a^+hrn,a + h.C.) 
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Fig. 1. Schematic diagrams of (a) the model Hamiltonian, (b) the 
band structure in the noninteracting case, and (c) a correspond- 
ing two-dimensional Fermi surface related to our ID model. 



(1) 



where c\ stands for the creation operator of an elec- 
tron with a spin a (=T, I) and an orbital m (= u,l) at 
site i and rii >mia — c| a Ci >mt<T . Here, t represents the 
hopping integral between the same orbitals and we set 
t = 1 in this study. The interaction parameters U, U' , 
J, and J' stand for the intra- and inter-orbital direct 
Coulomb interactions, exchange (Hund's rule) coupling, 
and pair-transfer, respectively. A denotes the energy dif- 
ference between the two atomic orbitals. For simplicity, 
we impose the relation J = J' . 

The model in eq. (1) is schematically shown in Fig. 
1(a). In a noninteracting case (U = U' = J = 0), 
the Hamiltonian eq. (1) yields dispersion relations rep- 
resenting the upper and lower band energies: e u (k) = 
2icos(fc) + y and e l (k) = 2tcos(fc) — — , where k is the 
wave vector. This band structure is schematically shown 
in Fig. 1(b). When the Fermi level Ek F crosses both 
the upper and lower bands, the system is metallic with 
electron and hole Fermi pockets corresponding to a char- 
acteristic band structure of the FeAs plane in iron oxyp- 
nictides, as shown in Fig. 1(c). 

We numerically diagonalize the model Hamiltonian up 
to 6 sites (12 orbitals) and estimate the Luttinger liquid 
parameter K p from the ground-state energy of finite-size 
systems using the standard Lanczos algorithm. 14 To re- 
duce the finite size effect, we impose the boundary con- 
dition (periodic or antiperiodic) on upper and lower or- 
bitals independently and chose both boundary conditions 



Fig. 2. K p as a function of U'(= 4J) for n = 5/3 (10 electrons/6 
sites) at A = 0.8, 1.6, 1.9, 2.3, and 2.8. The singlet ground state 
changes into a partially polarized ferromagnetic (S=l) state at 
U' ~ 2.5, 3.2, and 4.1 for A = 1.9, 2.3, and 2.8, respectively. The 
inset shows the energy difference Eo(<j>) — Eq(0) as a function of 
the external flux <j> for n = 2/3 (6 electrons/9 sites) at A = 1.2. 



to minimize \K p — 1|, where K® represents the K p of the 
finite-size system in a noninteracting case. The typical 
deviation of K p from unity becomes about ~ 0.1 for a 
six-site system. For simplicity, we will redefine K p as a 
renormalize value calculated using K p /K p , hereafter. 

On the basis of the Tomonaga-Luttinger liquid the- 
ory, 23-27 various types of correlation functions are deter- 
mined by the single parameter K p in the model which 
is isotropic in spin space. For a single-band model with 
two Fermi points, ikp, the SC correlation function de- 

cays as ~ r k p , while the CDW and SDW corre- 
lation functions decay as ~ cos(2fci?r)r~( 1+Ap '. Thus, 
the SC correlation is dominant for K p > 1, while the 
CDW or SDW correlation is dominant for K p < 1. On 
the other hand, for a two-band model with four Fermi 
points, i.e., ±kp 1 and ikp 2 , low-energy excitations are 
given by a single gapless charge mode with a gapped 
spin mode. 26-28 In this case, the SC and CDW correla- 
tions decay as ~ r 21< p and ~ cos[2(fci? 2 — kp 1 )r]r~ 2Kp , 
respectively, while the SDW correlation decays expo- 
nentially. Hence, the SC correlation is dominant for 
K p > 0.5, while the CDW correlation is dominant for 
K p < 0.5. In either case, the SC correlation increases 
with the exponent K p , and then K p is regarded as a 
good indicator of superconductivity. 29 As the noninter- 
acting K p is always unity, we assume that the condition 
of Kp > 1 for our model corresponds to the supercon- 
ducting state realized in oxypnictide superconductors. 

3. Phase Diagram 

Figure 2 shows K p as a function of U' for several val- 
ues of A at an electron density n — 5/3 (10 electrons/6 
sites), where we set J = U'/A with U = U' + 2 J. When 
U' increases, K p decreases for a small U', while it in- 
creases for a large U' in the case of A > 1.9, and then 
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Fig. 3. Phase diagram of the ground state with K p on the U' — J 
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Fig. 4. Phase diagram of the ground state with K p on the U — V 
parameter plane with J = U'/A for n = 5/3 (10 electrons/6 
sites) at A = 1.9. 



becomes larger than unity for U' > 2.3 in the case of 
A = 1.9. When J(= U'/A) is larger than a certain crit- 
ical value, the ground state changes into the partially 
polarized ferromagnetic state with a total spin S = 1 
from the singlet state with S = 0. We find that su- 
perconductivity is most enhanced in the vicinity of the 
partially polarized ferromagnetic state. To confirm su- 
perconductivity, we calculate the energy difference of the 
ground state, E {<p) — E o (0), as a function of the external 
flux (/>. As shown in the inset of Fig. 2, anomalous flux 
quantization is clearly observed for A = 1.9 while not for 
A = 0.8. 

In Fig. 3, we show the phase diagram of the ground 
state on the U' — J parameter plane under the condition 
of U = U' + 2J for n = 5/3 (10 clcctrons/6 sites) at 
A = 1.9. It contains the singlet state with S — to- 
gether with partially polarized ferromagnetic states with 
5 = 1 and 5 = 2. The singlet state with K p > 1, where 
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Fig. 5. Schematic diagrams of various types of superconducting 
pairing symmetries, S], n , 5" n , 5„ n , S^ n , and S^" with spin sin- 
glet pairings and T„ n , T" n , and with spin triplet pairings. 



we call it the SC phase, appears near the partially po- 
larized ferromagnetic region at J Z U'. It extends from 
the attractive region (V < 0) to the realistic parameter 
region with J ~ U'/A > 0, which is expected to corre- 
spond to that in the case of iron oxypnictides. 8 We have 
confirmed that similar phase diagrams are also obtained 
for A = 2.3 and 2.6. 

Figure 4 shows the phase diagram of the ground state 
on the U — U' plane under the condition of J = U'/A 
for n = 5/3 (10 electrons/6 sites) at A = 1.9. We ob- 
serve two types of SC phases with K p > 1, the SC I for 
U > U' and the SC II for U < U', in the wide param- 
eter region including both weak- and strong-correlation 
regimes. Note that the SC I corresponds to the SC phase 
shown in Fig. 3 and belongs to the realistic parameter 
region mentioned before. 

4. Pairing Correlation 

To examine the nature of these SC phases, we calculate 
SC pairing correlation functions for the various types 
of pairing symmetries schematically shown in Fig. 5. 
Explicit forms of the SC pairing correlation functions 
C(r) are given by 



S L( r ) = jjYl( c lu c lu Ci + r > l >l- Ci + r > l >^' 

i 

^on( r ) = ^ y^A C l.uA C i,u.l C i+r,u,l C i+r,u^), 
i 

S\mi r ) = ply X^( c l',T c !+i,U ~ c !u c !+i.z.t) 

i 

S"(r) = — V((ct ,4, ,-ct .ct,, t ) 

i 

X (ci+r+l,ti,|.C-i4-r,u,"|" Cz+r+l,u,'|'Ci+r, , u,.|.))5 
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where ^ n W, S U », S' n (r), S» n (r), and 5^ n (r) denote 
the singlet pairing correlation functions on the same site 
in the lower orbital, on the same site in the upper or- 
bital, between nearest-neighbor sites in the lower or- 
bital, between nearest-neighbor sites in the upper or- 
bital, and between the lower and upper orbitals on the 
same site, respectively. Furthermore, T^ n (r), T£ n (r), and 
T^(r) are the triplet pairing correlation functions be- 
tween nearest-neighbor sites in the lower orbital, between 
nearest-neighbor sites in the upper orbital, and between 
the lower and upper orbitals on the same site, respec- 
tively. 

In Fig. 6, we show the absolute values of various types 
of SC pairing correlation functions C(r)| for n = 5/3 (10 
electrons/6 sites) at A = 1.9, U' = 4J = 1.0, and U = 
—0.4. Here, the electronic state of the system belongs to 
the SC II phase, although the phase diagram for U < 
is not explicitly shown in Fig. 4. Note that \T^ n (r)\ < 
10~ 4 and S^(r = 3) = T^(r = 3) = 0, which arc not 
shown in Fig. 6. We find that S^ n (r) and 5" n (r) decay 
very slowly as functions of r, and |iS£ n (r = 3)| is the 
largest among the various \C(r = 3)| values. Therefore, 
a relevant pairing symmetry for the SC II phase seems 
to be the spin singlet pairing in the upper orbital band, 
mainly consist of 'on-site' pairing. It is considered that 
such a pairing in attractive region with U < is due 
to the intra-orbital attraction U . On the other hand, 
in repulsive region with U' > U > 0, the pairing may 
be due to charge fluctuation, which is enhanced by a 
large inter-orbital repulsion U' similarly to that in the 
case of the d-p model in the presence of the inter-orbital 
repulsion U p d- 30 

Next, we discuss the superconductivity in the SC I 
phase including the realistic parameter region mentioned 
before. Figure. 7 shows the absolute values of vari- 
ous types of SC pairing correlation functions \C(r)\ for 
n = 5/3 (10 electrons/6 sites) at A = 1.9, U' = 4 J = 1.0, 
and U = 2.4, where the system belongs to the SC I phase, 
as shown in Fig. 4. Here, |T n u n (r)|, \S l £(r = 3)|, and 

Ton(Y = 3)| are not shown, because these correlation 
functions are very small or zero. We find that |5" n (r)| is 
considerably suppressed compared with |5" n (r)| in con- 
trast to that in the case of the SC II phase. Furthermore, 

S^r) increases with increasing r except at r = 2. 




Fig. 6. Absolute values of various types of SC pairing correlation 
functions |C(r)| as functions of r for n=5/3 (10 clcctrons/6 sites) 
at A = 1.9, U'(= 4J) = 1.0, and U = -0.4, corresponding to 
the SC II phase. 



Therefore, the relevant pairing symmetry for the SC I 
phase seems to be an extended spin singlet pairing, and 
mainly consist of nearest-neighbor site pairing. 




Fig. 7. Absolute values of various types of SC pairing correlation 
functions |C(r)| as functions of r for n=5/3 (10 clcctrons/6 sites) 
at A = 1.9, U {= 4J) = 1.0, and U = 2.4, corresponding to the 
SC I phase. 

Recently, weak-coupling approaches such as RPA 
and perturbation expansions have shown that the sign- 
reversing s-wave (s±-wave) pairing is realized in iron 
oxypnictide superconductors. 9-13 The order parameter 
of such a pairing is considered to change its sign be- 
tween hole and the electron Fermi pockets. To compare 
our result with the result obtained by weak-coupling ap- 
proaches, we examine the SC pairing correlation function 
between the lower and upper orbitals, such as 

Sin U (r) = 2^ E < A U^Kn(i + r)> 
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with 

A™ n {tf = c I,m,T C l+l,m,| ~ c I,ro,J. C I+l,ro,T ( m = ^ M )' 

We also define T^~ u (r) as well as S„~ u (r) in the above 
equation. 

When the s±-wave pairing is dominant, the values of 
the interorbital SC pairing correlation function are ex- 
pected to be negative, since the Fermi surface of the 
lower (upper) orbital band in our model corresponds to 
a hole (electron) Fermi pocket, as shown in Fig. 1(b). In 
Fig. 8, we show the interorbital pairing correlation func- 
tions ^""(r) and T^~ u (r)(see also inset) for the same 
parameters in Fig. 4 corresponding to the SC I phase. 
We see that the values of T„~ u (r) are positive and very 
small, while those of 5„~ u (r) are negative except at r = 3 
and not so small. This result suggests that the rele- 
vant pairing symmetry of the SC I phase is the spin- 
singlet s±-wave pairing, which agrees with the result of 
the weak coupling approaches. Therefore, we expect that 
the s±-wave pairing proposed on the basis of the result of 
weak-coupling approaches is realized in the wide param- 
eter region including both weak- and strong-correlation 
regimes. 



U=2.4 U'=4J=1.0 




I I i I i I i 

1 r 2 3 

Fig. 8. Pairing correlation functions Snn"^) and T n n U (r). Here, 
we show the absolute value of the correlation functions at U = 
2.4, A = 1.9, and U'(= 4J) = 1.0 for n=5/3 (10 electrons/6 
sites). The inset shows a schematic diagram of the pairing sym- 
metries, Sin U (»*) an d Tnn U (»*)- 

Finally, we discuss the mechanism of superconductiv- 
ity in the SC I phase. The s±-wave pairing is consid- 
ered to be mediated by the antifcrromagnctic fluctuation 
due to the nesting effect between electron and hole pock- 
ets. 9-13 At first glance, the SC I phase is located adjacent 
to the partial ferromagnetic phase (S=l), and then, the 
ferromagnetic fluctuation seems to be related to super- 
conductivity. To examine the relationship between spin 
fluctuation and superconductivity, we also calculate the 
spin correlation function for finite-size systems, where a 
short-range spin correlation is considered to be crucial 
for the superconductivity. We obtain the ferromagnetic 
and antiferromagnetic components of the spin correla- 
tion as a function of U'(= 4 J) for a fixed U — 1.5 (see 
also Fig. 4) and find that antiferromagnetic (ferromag- 



netic) correlation increases (decreases) with decreasing 
U' together with increasing K p (not shown). Therefore, 
we conclude that the antiferromagnetic spin fluctuation 
is responsible for the s±-wave pairing in the SC I phase. 

5. Summary and Discussion 

We have investigated the superconductivity of the one- 
dimensional two-orbital Hubbard model in the case of 
electron and hole Fermi pockets corresponding to a char- 
acteristic band structure of iron oxypnictide supercon- 
ductors. To obtain reliable results including those in the 
strong-correlation regime, we used the exact diagonal- 
ization method and calculated the critical exponent K p 
on the basis of the Luttinger liquid theory. It has been 
found that the system shows two types of SC phases, 
the SC I for U > U' and the SC II for U < U', in the 
wide parameter region including both weak- and strong- 
correlation regimes. 

We have also calculated various types of SC pairing 
correlation functions in the realistic parameter region of 
the iron oxypnictides. The result indicates that the most 
dominant pairing for the SC I phase is the intersite in- 
traorbital spin-singlet with sign reversal of the order pa- 
rameters between two Fermi pockets. The result is con- 
sistent with the sign-reversing s±-wave pairing that has 
recently been proposed on the basis of the result obtained 
by weak-coupling approaches for iron oxypnictide super- 
conductors. This indicates that the s±-wave pairing is 
realized not only in a weak-correlation regime but also 
in a strong-correlation regime. We have also calculated 
the spin correlation function and found that antiferro- 
magnetic spin fluctuation is responsible for the s±-wave 
pairing in the SC I phase. 

As for the SC II phase, the most dominant pairing 
is found to be the on-site intraorbital spin-singlet pair- 
ing, which is consistent with the ordinary s-wave pairing 
of BCS superconductors. However, the superconducting 
mechanism of this phase is due to the charge fluctuation 
enhanced by the interorbital Coulomb interaction and is 
different from the conventional BCS superconductivity 
due to the electron-phonon interaction. Although the 
SC II phase seems to be realized only for the unrealis- 
tic parameter region in our model, it might be realized 
for a realistic parameter region in the d-p model, which 
is closer to iron oxypnictides. 13,30 We will address such 
a problem by applying the present method in the d-p 
model in the future. 
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